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Abstract. We present a detailed analysis of time-resolved optical spectra of the ZZ Ceti white dwarf, HS 
0507-I-0434B. Using the wavelength dependence of observed mode amplitudes, we deduce the spherical degree, 
£, of the modes, most of which have i — 1. The presence of a large number of combination frequencies (linear 
sums or differences of the real modes) enabled us not only to test theoretical predictions but also to indirectly 
infer spherical and azimuthal degrees of real modes that had no observed splittings. In addition to the above, 
we measure line-of-sight velocities from our spectra. We find only marginal evidence for periodic modulation 
associated with the pulsation modes: at the frequency of the strongest mode in the lightcurve, we measure an 
amplitude of 2.6 ± l.Okms"^, which has a probability of 2% of being due to chance; for the other modes, we find 
lower values. Our velocity amplitudes and upper limits are smaller by a factor of two compared to the amplitudes 
found in ZZ Psc. We find that this is consistent with expectations based on the position of HS 0507-I-0434B in the 
instability strip. Combining all the available information from data such as ours is a first step towards constraining 
atmospheric properties in a convectionally unstable environment from an observational perspective. 

Key words, stars: individual: HS 0507+0434, ZZ Psc - white dwarfs - oscillations - convection 



1. Introduction 

The apparent simplicity of white dwarfs belies more com- 
plex and ill-understood processes occurring in their inte- 
riors. Pulsations, however, offer the precious possibility of 
probing the interiors of these stars thereby yielding not 
only fundamental stellar parameters but also important 
clues of their prior history which can, in principle, be re- 
constructed from their present interior structure. 

Along the white dwarf cooling track, there are three 
regions of instability: at T^s ^ 110 kK, T^s ^ 24kK, and 
Toff - 12 kK, populated by the GW Vir, V 777 Her (DBV), 
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and ZZ Ceti (DAV) types, displaying strong lines in their 
optical spectra of Heii and Civ, Hei, and H respectively. 

Given the relatively long pulsation periods of white 
dwarfs, the realisation that white dwarfs are non-radial 
gravity-mode pulsators (Chanmugam 1972; Warner & 
Robinson 1972) came soon after the discovery of the first 
variable white dwarf (Landolt 1968); that their photomet- 
ric variations are primarily a manifestation of temperature 
perturbations rather than due to variations in geometry 
or surface gravity came only a decade later (Robinson et 
al 1982). 

In order to subject the ZZ Cetis to asteroseismologi- 
cal analysis and to provide constraints for pulsation mod- 
els, it is crucial that the eigenmodes associated with the 
observed periodicities be identified. Observationally, this 
means determining the value of the spherical degree, i, 
and azimuthal order, m. A third quantity, not an observ- 
able, is the radial order n; it specifies the number of nodes 
in the radial direction and can only be inferred by detailed 
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comparison of observed mode periods with those predicted 
by pulsation models. 

Mode identification for the ZZ Cetis is fraught with 
difficulties. In part, the analysis has been hampered by 
an insufficient number of pulsation modes excited to ob- 
servable levels and mode variability over several different 
time scales. For the vast majority of ZZ Cetis, results have 
remained somewhat ambiguous as the prerequisites for as- 
tcroseismological analysis were not met. Even the Whole 
Earth Telescope (WET) campaigns (Nather et al. 1990) 
on several objects (e.g. G 117-B15A, ZZ Psc, Kepler et al. 
1995; Kleinman et al. 1998) were thwarted cither by the 
small number of modes and lack of clear multiplet struc- 
ture exhibited by these objects. Most efforts have focused 
on identifying similarities between the pulsational spectra 
of different stars to constrain the mode identification. This 
ultimately has a bearing not only on the determination of 
fundamental stellar properties but also on the mass of the 
superficial hydrogen layer (e.g. Bradley 1998). It is clear 
that there is an acute need for more direct and comple- 
mentary methods of pinning down the identification of the 
eigenmodes. 

To this end, Robinson et al. (1995) presented a method 
for inferring ^ based on the the wavelength-dependence 
of limb darkening, due to which observed mode ampli- 
tudes vary with wavelength in a manner that depends on 
t, but not on any of the other properties of the pulsation 
mode, such as m or amplitudc"'^ . Thus, in a given star, 
modes having the same spherical degree will behave in 
the same manner. Robinson et al. (1995) acquired photo- 
metric data in the ultraviolet of the ZZ Ceti star G 117- 
B15A. Application of Bayes' theorem and quantitative use 
of model wavelength-dependent pulsation amplitudes led 
them to infer that the largest amplitude mode of G 117- 
B15A had ^ = 1. 

More recently, van Kerkwijk et al. (2000) and Clemens 
et al. (2000) used a variant of the above method to identify 
the spherical degree of the pulsation modes observed in ZZ 
Psc (a.k.a. G 29-38) - a star that has been notoriously er- 
ratic in the pulsation modes that it excites. Their investi- 
gation, which was based on amplitude changes within the 
Balmer lines at visual wavelengths only, yielded empirical 
differences between the modes that were best interpreted 
as several £ = 1 modes and one £ = 2 mode; the presence 
of modes of differing £ obviated the need for quantitative 
model comparisons making their identification more se- 
cure than any previous attempt. A surprising by-product 
of their analysis was the detection of variations in the line- 
of-sight velocity associated with the pulsations. Given the 
instrumentation available at that time, these were thought 
to be too small to measure (Robinson et al. 1982). 

The measurement of line-of-sight velocity variations 
associated with the pulsations can be used in two com- 
plementary ways: (i) they provide a means with which to 
verify and constrain the theories of mode driving and, (ii) 

^ At least for pulsation amplitudes of up to ~5% (Ising & 
Koester 2001). 



under certain theoretical assumptions, velocity variations 
in conjunction with flux changes provide an important 
new tool with which to probe the outer layers of pulsating 
white dwarfs. 

In this study, we will only attempt to interpret our 
observations within the context of theories of mode driv- 
ing via convection (Brickhill 1983, 1991; Goldreich & Wu 
1999a,b), as opposed to theories which purport mode 
driving by variants of the classical K-mcchanism (e.g. 
Dzicmbowski & Koester 1981; Dolez & Vauclair 1981; 
Winget et al. 1982). Unfortunately, as far as we are aware, 
testable predictions for these models are not, as yet, avail- 
able. 

Within the context of the "convective-driving" pic- 
ture the convection zone responds to perturbations from 
the adiabatic interior on time scales very much shorter 
(~ Is) than the mode periods (hundreds of seconds). 
This allows the convection zone to absorb and release 
the flux perturbations cyclically, thus driving the pulsa- 
tions. Combination frequencies (linear sums and differ- 
ences of the real modes) arise naturally in the above pic- 
ture and can, in principle, provide additional information. 
Furthermore, the relation between the flux and velocity 
variations allows one to estimate the total thermal capac- 
ity of the convection zone. Qualitative relations derived 
by Goldreich & Wu (1999a,b) show that these are also 
sensitive to £. We explore these issues in Sect. 7. 

Based on the success of identifying the £ index of 
the modes of ZZ Psc using time-resolved spectroscopy, 
we present here, results from a similar analysis on HS 
0507-I-0434B. Our primary aims are to attempt to de- 
termine the spherical degree (£) of the pulsations and to 
search for line-of-sight velocity variations associated with 
the pulsations. We also hope to constrain the properties 
of convection, a process which is poorly understood and 
therefore one of the main sources of uncertainty in the 
models. 

We begin with a brief introduction to HS 0507-1-0434 

in Sect. 2 followed by a description of the data in Sect. 3. 
In Sects. 4, 5, and 6, we attempt to extract the flux and 
velocity amplitudes and phases and the change in pulsa- 
tion amplitudes with wavelength. From Sect. 7 onwards, 
we apply the constraints provided by all of the above to 
extricate quantities that allow us to both test and subse- 
quently use the theory. 

2. HS 0507-1-0434 

HS 0507-1-0434, like many other white dwarfs, is the by- 
product of a survey for faint blue objects (the Hamburg 
Quasar Survey in this case). It comprises two DA white 
dwarfs in a common proper-motion pair which, as dis- 
cussed by Jordan et al. (1998), is of particular interest for 
two reasons. First, the effective temperatures of 20 kK for 
HS 0507A and 11.7kK for HS 0507B are such that the at- 
mosphere of the former is completely radiative, while the 
latter has an outer convection zone. Thus, one can hope to 
use the parameters inferred from HS 0507A, which should 
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Fig. 1. Sample and average spectra of HS 0507A and B in 
the 3650-5800 A range with the Balmer hnes indicated. The 
average spectra are offset from the sample spectra by f mjy 
and 0.65 mJy for the A and B components respectively. 



be secure - since radiative atmospheres consisting of pure 
hydrogen are well understood - to calibrate models for HS 
0507B, which have to rely on uncertain models for convec- 
tion. Jordan et al. (1998) found that only with relatively 
inefficient versions of the Mixing Length Theory (MLT) 
were they able to find consistent solutions for both com- 
ponents in the system. 

The second reason that HS 0507+0434 is of interest 
is that the temperature of HS 0507B as derived from 
its spectrum - placed it squarely within the ZZ Ceti in- 
stability strip; fast photometry on this object revealed it, 
as expected, to be variable (Jordan et al. 1998). A study 
of the temporal behaviour of HS 0507-I-0434B, based on 
photometric measurements collected over a total of seven 
consecutive nights, was recently carried out by Handler 
& Romero-Colmenero (2000). These authors were able to 
resolve a set of three equally spaced triplets^. Under the 
assumption that the splitting was due to the effect of slow 
rotation on £ = 1 modes, the observed triplets yielded an 
estimate of the rotation rate (~ 1.5 days) and the m val- 
ues of the multiplets. They also found that in all three 
triplets, the m = component was much weaker than the 
m = — 1 and m = +1 components, which were of roughly 
similar strength. Assuming that the intrinsic amplitudes 
were the same for all m components, they estimated the 



Assuming spherical symmetry, frequencies of modes hav- 
ing the same n and £ are degenerate. Slow rotation and/or a 
magnetic field lifts this degeneracy resulting in 2£ + 1 and £ + 1 
split components respectively. 



inclination of the rotation axis with respect to the line of 
sight to be ~ 79°. The independent £ identifications af- 
forded by clearly split multiplets are highly desirable as 
they can help to validate the identifications that rely on 
time-resolved spectroscopy only, given the differences be- 
tween model spectra and observations. 

As will become clear in the following sections, the ad- 
vantage of having a flux and velocity reference in the same 
slit as the target greatly increases the accuracy of subse- 
quent measurements by making it possible to not only 
divide out atmospheric fluctuations, but also to ensure 
that small random movements of the target in the slit are 
accounted for in the determination of the Doppler shifts 
of the Balmer-lines, thus making HS 0507-1-0434 an ideal 
system on which to test theoretical predictions. 



3. Observations and Data reduction 

HS 0507+0434 was observed on the nights of 10*i^ and 
11*^ December 1997 using the Low Resolution Imaging 
Spectrometer (LRIS) on the Keck II telescope (Oke et 
al. 1995). An 8'.'7-wide slit was used together with a 
600 linemm"^ grating covering the 3450-5960 A range at 
1.25 Apixel"^. On-chip binning by a factor of two was 
implemented in the spatial direction. The first night was 
photometric and the seeing of 1'.'2 set the wavelength res- 
olution to 7 A. There were thin patches of cirrus on the 
second night; furthermore, there was a great deal of wind- 
shake, which resulted in a worse resolution (9 A) although 
the situation improved towards the end of the run as the 
elevation decreased. A contiguous set of 560 24 s expo- 
sures were acquired from 7:13:35 to 13:17:17 U.T. on the 
first night and a further 280 exposures were obtained from 
7:10:50 to 10:12:24 U.T. on the following night. For both 
nights, the series was preceded by exposures of the flux 
standard G 191-B2B (five and ten integrations of 4 s re- 
spectively), followed by HgKr arc spectra for wavelength 
calibration, and halogen flat-field frames. A total of 64 
frames from the second night had to be discarded due to 
a malfunction of the guidcr, which resulted in the tar- 
gets drifting off the slit. We decided to fix the slit at the 
position angle required to have both components of the 
system in the slit and deal with the effects of differential 
refraction during the reduction of the data. The decision 
to use a wide slit stemmed from the need to acquire 'pho- 
tometric' measurements of HS 0507 A and B in addition 
to measurements of variations in the line-of-sight velocity 
(see Sect. 4). 

The reduction of the data was carried out using 
MIDAS^ and using routines running in the MIDAS envi- 
ronment written specifically for the LRIS instrument and 
this data set. In detail, our reduction procedure entailed: 
(i) subtracting the CCD bias level of each amplifier sep- 
arately using the overscan region; (ii) correcting for the 
(small) uou-liiiearities that most CCDs are prone to, using 



The Munich Image Data Analysis System, developed and 
maintained by the European Southern Observatory. 
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Fig. 2. Light and velocity curves of HS 0507B, shown here for 
the first night only for the sake of clarity, t' = t— 10:15:14 U.T. 
i.e. relative to the midpoint of the time series, (a) Fractional 
variations in the continuum region between ~ 5300 - 5800 A. 
The variations are with respect to the continuum of HS 0507A 
in the same wavelength region. Note that the maxima are 
stronger and sharper than the minima as is typical for ZZ Ceti 
type pulsators. The bottom curve shows the residuals (offset 
by +0.92) after fitting 26 sirmsoids with amplitudes larger than 
0.15% (see Table 1); (b) variations in the line-of-sight velocity 
also with respect to HS 0507A - as derived from fitting pro- 
files to H/3, H7, and H(5. No obvious variations are present and 
indeed we detect only marginal velocity variations associated 
with the pulsation modes (Sect. 5). 



linearity curves of the CCD derived from previous obser- 
vations using the LRIS instrument; (iii) correcting for the 
gain difference between the 2 amplifiers as derived from 
halogen lamp frames; (iv) flat-fielding using an average of 
the halogen frames for the first night and an average of 
dome-flats for the second night (as these gave smoother re- 
sults); (v) sky subtraction; (vi) correcting for the error in- 
troduced by dividing by a flat field taken through slightly 
non-parallel slits (described in more detail below); (vii) 
optimal extraction of the spectra using a method akin to 
that of Horne (1986); (viii) wavelength calibration using 
Hg/Kr frames (by tabulating wavelengths for each pixel 
rather than by rebinning the spectra), with an offset de- 
termined from star A, and including a correction for re- 
fraction (see below) ; and (ix) fiux calibration with respect 



to the fiux standard G 191-B2B, using the model fiuxes 
of Bohlin et al. (1995) and using the extinction curve of 
Beland et al. (1988) to correct for small differences in air- 
mass. 

Among the usual preprocessing stages described above, 
two additional non-standard steps were required given our 
use of a wide slit, namely the need to correct for non- 
parallel slit-jaws and the need to account for the effects 
of random stellar wander in the slit due, for instance, 
to sceing-related changes, windshake, and tracking uncer- 
tainties. 

The former correction is necessitated by the fact that 

while the flat-fields and sky background are influenced by 
the shape of the slit, the light from the targets is not, 
as no light falls outside the slit jaws. The sky background 
was determined by fltting flrst degree polynomials for each 
step along the dispersion direction excluding two regions 
of 26 and 30 binned pixels corresponding to HS 0507A 
and HS 0507B respectively, and multiplying the result by 
a function describing the shape of the slit. Finally, we di- 
vided the images by an average flat field frame that was 
normalised along the spatial direction in order to main- 
tain the correct relative fiux ratio between the two target 
objects. 

A second non-standard correction is required to purge 
the data of any extraneous effect that introduces shifts 
in the positions of the Balmer lines, such as instrumen- 
tal fiexure, refraction, and wander. In principle, all these 
changes should affect HS 0507A and B in the same way, 
and one could simply measure wavelengths relative to HS 
0507A only. However, we prefer the following methodol- 
ogy, as it not only permits us to assess the quality that 
one may expect in the absence of such a fortuitous local 
calibrator, but it also permits us to verify the stability of 
HS 0507A and provides insight into any systematic effects 
that may have crept into our analysis of other targets. 
We first correct the spectra of both stars for instrumental 
fiexure and refraction'* and only then apply a correction 
for stellar wander, using the measured offsets of H7 of HS 
0507A as fiducial points. Having thus corrected each pixel 
in the spectrum for the effects of instrumental fiexure, 
differential refraction, and random wander in the slit, we 
redid the wavelength and fiux calibrations. Any remaining 
offsets of the Balmer lines from their respective laboratory 
wavelengths must now be due to intrinsic variations only. 

An estimate of the remaining scatter due to wander in 
the derived velocities can be obtained from the shifts in 
the slit of the spatial profiles of both HS 0507A and HS 
0507B and assuming that the scatter in these shifts is also 
representative of the scatter in the dispersion direction. 
Having fitted Gaussian functions to the spatial profiles, 
we take the standard deviation of the difference of the 
resulting centroid positions of HS 0507B with respect to 
those of HS 0507A as an indicator of the scatter and obtain 
a vi\h\v of ~2kin s~^at H7 as an average value for both 



^ See Appendix B where we scrutinize HS 0507A for any 
signs of variability. 
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Fig. 3. (a) Fourier transform of the Hght curve (top) and residuals (bottom) offset by —0.2%. The maximum amplitude of 
peaks longward of 9000/iHz is ~ 0.09%. (b) Fourier transform of the velocity curve. The frequencies inferred from the light 
curve are indicated. None of the frequencies present in the light curve have obviously significant associated velocity modulation. 
This is confirmed by our detailed analysis (Sect. 5), with the possible exception of F2, for which we find that associated signal 
has a probability of only 2% of being due to chance, constituting a marginal detection. (Note that without prior information 
from the light curve regarding which frequencies to expect, this peak would also be insignificant). Both (a) and (b) have been 
constructed using data from both nights. 



nights. This is well below the measurement errors (see 
Sect. 5). Note that the root-mean-square scatter for ZZ 
Jr^sc was 14kms-iat H7 (van Kerkwijk at al. 2000). The 
value of having a reference object in the slit thus cannot 
be overemphasizeid. 

The final average spectra of HS 0507+0434 A and B, 
as well as sample individual spectra, are shown in Fig. 1. 

4. Periodicities in the Light Curve 

Continuum light curves for each night were constructed 

by dividing the average flux in the 5300-5800 A region of 
HS 0507B by that of HS 0507A (see Figs. 1, 2a). Note 



that the time axis has been computed relative to the mid- 
dle of the time series in order to minimize any covariance 
between the frequencies and phases of the sinusoids we 
subsequently use to fit the light curve {t' = t—16 : 33 : 20 
U.T.). The Fourier transform of the light curve obtained 
by merging the data from both nights is shown in Fig. 3a 
and is typical of that of the ZZ Cetis. 

We determined the periodicities consecutively in terms 
of decreasing amplitude by means of an iterative process: 
using an approximate value for the frequency of the high- 
est peak, we fitted a function of the form A cos(27r/i' — 
4>) + C where A is the amplitude, / the frequency, </> the 
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Table 1. Pulsation frequencies and other derived quantities from the light and velocity curves of HS 0507B 



Mode Period Frequency 3>i Ay *v Rv A4>v 
(s) iflUz) n (%) (kms-^) (°) (Mmrad-^) (°) 



Real Modes: 



Fl (/lo) m = 1 


354.9 ± 0.03 


2817.5 ± 0.3 


-27 ± 4 


0.68 ± 0.11 


2.0 ± 1.0 


61 ± 29 


17 ± 


9 


88 ± 29 


F2 (/s) rn = -1 


355.8 ± 0.02 


2810.9 ± 0.2 


-57 ± 1 


2.27 ± 0.06 


2.6 ± 1.0 


-13 ± 23 


6 ± 


3 


44 ± 23 


F3 (/2) rn = -1 


557.7 ± 0.03 


1793.2 ± 0.1 


-6 ± 1 


1.87 ± 0.04 


0.8 ± 0.8 




4 ± 


4 




F4 (fl) 


743.0 ± 0.13 


1346.1 ± 0.2 


-153 ± 2 


1.39 ± 0.04 


0.3 ± 1.1 




2 ± 


9 




F5 (/a) m = -1 


446.2 ± 0.05 


2241.1 ± 0.2 


-0.2 ± 3 


1.10 ± 0.09 


1.6 ± 1.1 




10 ± 


7 




F6 (/t) m = 1 


444.8 ± 0.06 


2248.4 ± 0.3 


74 ± 3 


1.36 ± 0.07 


2.3 ± 1.1 


79 ± 26 


12 ± 


5 


4 it 27 


F7 


286.1 ± 0.04 


3495.8 ± 0.5 


147 ± 6 


0.36 ± 0.04 


0.8 ± 0.8 




9 ± 


9 




itiozticizzoH' ivi odest 




















Mode 


Period 


Frequency 


3" I, 


Al 


Av 




Rc 








(s) 


(/iHz) 


( ) 


(%) 


(knis~^) 


( ) 






( ) 


F4+F6 


278.2 


3594.4 


— 131 ± 4 


U.oU ± U.Uo 


0.9 it 1.1 




20.7 lb 


1.8 


—53 ± 5 


F1+F3 


216.9 


4610.7 


—32 ± 7 


0.38 ± 0.04 


0.2 ± 1.0 




15.0 it 


3.0 


1 ± 8 


F2+F3 


217.2 


4604.2 


— 73 ± 8 


0.30 ± O.Oo 


1.1 it 1.0 




3.5 it 


0.6 


— 10 lb 8 


F2-|-F4 


Z4:U.O 


4156 8 


129 i 11 


fl 9n -1- n n/1 

U. ZU IE U . U4i 


4 i 8 




O.U HI 


u.o 


21 ^ 12 


F3+F6 


247.4 


4041.7 


71 ± 8 


0.33 ± 0.04 


1.0 ± 1.0 




6.4 ± 


0.9 


2 ± 8 


F3+F5 


247.9 


4034.4 


-57 ± 17 


0.16 ± 0.04 


0.4 ± 1.1 




3.9 it 


1.1 


-51 it 17 


F1+F5 [F2+F6] 


197.7 


5058.6 


-12 ± 4 


0.50 ± 0.04 


1.4 ± 0.8 




32.8 it 


6.2 


15 it 7 


F2-F4 


682.6 


1464.9 


41 ± 6 


0.37 ± 0.04 


1.3 ± 0.8 




5.9 it 


0.6 


-55 it 6 


F5+F6 


222.7 


4489.6 


80 ± 8 


0.27 ± 0.04 


0.8 ± 0.8 




9.0 it 


1.5 


5 it 9 


F1+F2 


177.7 


5628.4 


-140 ± 10 


0.21 ± 0.04 


0.5 ± 0.8 




6.7 it 


1.6 


-57 it 11 


F3-F4 [F5-F3] 


2236.5 


447.3 


-102 ± 9 


0.23 ± 0.04 


0.3 ± 0.8 




4.3 it 


0.7 


111 it 10 


F1-F5 


1735.0 


576.4 


144 ± 10 


0.22 ± 0.04 


0.9 ± 0.8 




14.4 it 


3.5 


170 ± 11 


F3+F4 


318.6 


3139.2 


160 ± 11 


0.19 ± 0.04 


1.0 ± 0.8 




3.6 it 


0.7 


-41 ± 11 


2F3 [F4+F5] 


278.8 


3586.5 


-128 ± 17 


0.19 ± 0.05 


1.4 ± 1.0 




5.4 it 


1.5 


-116 it 17 


F1+F2+F3 


134.7 


7421.6 


-128 ± 14 


0.16 ± 0.04 


0.3 ± 0.8 




90 it 


30 


-38 it 14 


F2+F4+F6 [2F3+F1] 


156.1 


6405.2 


151 ± 7 


0.28 ± 0.04 


0.7 ± 0.8 




108 it 


16 


-74 ± 8 


F3+F4+F6 


185.6 


5387.6 


152 ± 8 


0.25 ± 0.04 


0.7 ± 0.8 




117 it 


19 


-124 it 9 


2F4+F6 


202.4 


4940.3 


-14 ± 13 


0.17 ± 0.04 


1.8 ± 0.8 




208 it 


50 


-143 it 14 


F4+F6-F5 [m=l, F4] 


739.0 


1353.3 


1 ± 18 


0.24 ± 0.05 


1.3 ± 1.0 




190 ± 


43 


79 ± 18 



Notes. The frequencies listed in brackets next to the real modes refer to the nomenclature preferred by Handler &z Romcro-Colmencro (2000). The 
curly brackets next to two pairs of real modes indicate that these are multiplets sharing the same i but differing m values (see also Sect. 7.2). 
We have followed the convention that positive values of m correspond to prograde modes. Rv — Ay /{27rfAL), in units of Mm rad~^, has the 
physically intuitive meaning of being the typical distance a fluid element travels on the surface for a pulsation with a fractional amplitude of unity. 

and are the phases of the sinusoids fit to the light and velocity curves respectively with A^>v = — ^l- For the combination frequencies, 
Rc = A^^/{nijAiA{) where nij is unity for i = j, 2 for 2-niode combinations etc. and 3 for 3-mode combinations involving the first harmonic, 
e.g. 2F4+F6. A^c = ^^l^^ ~ (^l =t 'J*;^). Potentially misidentified modes are listed in square brackets next to the combination modes. Note that 
this list is not an exhaustive one. The frequency of the combination F4+F6— F5 is at the location of the m = 1 multiplet c.f. F4 (we infer in Sect. 
6 that F4 has ^ = 1 and in Sect. 7.2 that it probably has m = —1). 



phase and C, a constant offset. The fit yielded A, f and 
(p. The process was repeated, adding a new sinusoid each 
time, and fitting for all parameters simultaneously until 
no peaks with amplitudes ^ 0.15% could be identified in 
the Fourier transform of the residuals. 

Combination frequencies were identified by searching 
for linear combinations of the real modes and these were 
fitted by fixing the frequencies to those of the correspond- 
ing combination of the real mode frequencies. We imposed 
the requirement that the combination frequency in ques- 
tion have the smallest amplitude of the frequencies in- 
volved. Identification of combination frequencies was not 
always easy due to degenerate combinations (e.g. F3— F4 
~ F5— F3 given our resolution). In such cases, we opted 
for the largest amplitude of the combination frequency 
and the lowest possible error yielded by the fit. The light 
curve is not free of periodicities after having subtracted 



the frequencies listed in Table 1; although we continued 
to fit further combination modes, we found that the choice 
of combinations became rather arbitrary and hence we do 

not report these here. 

Two statements can already be made on the basis of 
Table 1. First, as Fl and F2 have the same splitting as 
the F5 and F6 multiplet components, they therefore very 
probably also share the same i value. From the analysis of 
Handler & Romero-Colmenero (2000), who detect triplets 
in each of these groups, they are likely to have £ ~ 1 
and TO = ±1. Second, the periods of the real modes are 
very similar to those observed by Kleinman et al. (1998) 
in a time series spanning 10 years for ZZ Psc, a similar, 
but slightly cooler white dwarf. In that star, these authors 
find real modes at periods ranging from 110 - 915 s. For 
£ = I, these correspond to successive radial orders (n) 
from 1 to 18. The modes in common with HS 0507B are 
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at 284, 355, and 730 s, which they identify as having ra- 
dial orders of 4, 5, and 13 respectively. The n = 7 mode 
in their model is the only mode (from n = 1 - 18) not 
detected in the ZZ Psc data. Interestingly, a mode at the 
expected frequency is present in HS 0507B as the F5,F6 
multiplet pair (445 s) at the expected separation for ^ = 1 
modes. It is tempting to conclude that, as for the hot ZZ 
Cetis (Clemens 1993), the cool ZZ Cetis have remarkably 
similar structure. One should bear in mind, however, that 
similar period structures may well be produced for a num- 
ber of different combinations of mass and hydrogen-layer 
thickness (e.g. Bradley 1998). 

We refrain from a detailed comparison with the study 
of Handler & Romero-Colmenero (2000) given our lower 
frequency resolution, but a few points are noteworthy. 
While the frequencies of the modes detected in our data 
are entirely consistent with those found by Handler & 
Romero-Colmenero (2000), our lower resolution docs not 
permit us to confirm the presence of the m = multi- 
plet components independently (see Fig. 4a); these were 
also the weakest modes in all triplets and thus should not 
greatly influence the amplitudes of the m = ±1 compo- 
nents. As a check, we attempted to recover the real modes 
identified by Handler & Romero-Colmenero (2000) that 
were not detected by us, by first imposing the amplitudes 
and frequencies listed by Handler & Romero-Colmenero 
(2000) and subsequently leaving the amplitudes free to 
vary. The residuals are shown in Fig. 4b. We find that our 
data are consistent with all but one of the components and 
have amplitudes consistent with those derived by Handler 
& Romero-Colmenero (2000). The only exception is the 
m = 1 mode at 1800.7 ^Hz (in the same multiplet as our 
F3), for which Handler & Romero-Colmenero (2000) found 
an amplitude of 1.7%, while it has an amplitude of only 
^ 0.5% in our data and is not even required in the fit. 
For completeness, we note that the amplitudes and phases 
we derive using the frequencies of Handler & Romero- 
Colmenero (2000) are entirely consistent with those listed 
in Table 1. 

5. Periodicities in the line-of-sight velocities 

Variations in the equivalent widths and profiles of the 
Balmer lines reflect temperature and line-of-sight veloc- 
ity variations. The details of the complex interactions be- 
tween flux and velocity that result in the observed spectral 
line are not well-understood, but are likely to be separable 
to flrst order. 

In order to search for line-of-sight velocity variations, 
we determined Doppler shifts for the Balmer lines by flt- 
ting a combination of a Gaussian and a Lorcntzian profile 
which provided a good fit to the lines; the continuum was 
approximated by a line. We chose the wavelength intervals 
4660-5088 A, 4214-4514 A, and 4028-4192 A for H/3, H7, 
and H(5 respectively to carry out the fits and insisted that 
the central wavelengths of the Gaussian and Lorentzian 
functions be the same. Wc repeated the above process for 
H7 only, fitting the logarithm of the fiux instead as a check 




I I I I I I I I I I I I I 

2000 2500 3000 

Frequency (//Hz) 

Fig. 4. (a) The top curve shows the FT of the Hght curve in 
the 1600 - 3000 /xHz range. The middle curve shows the FT 

of the residuals after subtracting sinusoids with frequencies 
corresponding to all real modes except Fl and F6, which are 
multiplets in the same period group as F2 and F5 respectively, 
while the bottom curve shows the effect of including Fl and F6 
in the fit. These two lower curves have been offset by —0.25% 
and —0.40% respectively. While there are some feeble signs of 
residual modulation around Fl and F2, this is clearly not the 
case for F5 and F6. Fitting an additional sinusoid near the 
F1,F2 pair in the hope of recovering the missing component is 
of no avail, (b) Top curve: residuals in the FT after subtracting 
three sets of triplets with frequencies and amplitudes fixed to 
those listed in Handler & Romero-Colmenero (2000). The large 
residual corresponds to the m = 1 mode at 1800.7/LtHz showing 
that the amplitude of this mode has changed considerably since 
our observations, from 1.7 to < 0.5%. The bottom curve (offset 
by —0.3%) shows the effect of leaving the amplitudes free to 
vary while leaving the frequencies fixed. 



of our (additive) fitting method. The velocities thus de- 
rived were nearly identical with those obtained by fitting 
the fiux only. Possible systematic effects arising from dif- 
ferent methods of deriving these velocities are discussed in 
detail in van Kcrkwijk ct al. (2000). The resulting veloc- 
ity curve is shown in Fig. 2b, while the associated Fourier 
transform is shown in Fig. 3b. The typical uncertainty in 
a single measurement for HS 0507B is 14kms~^, much 
larger than those due to differential wander between HS 
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0507A and B which was estimated to be ^ 2kms ^(Sect. 
3). 

Simple checks confirm that the errors in the measure- 
ment of the linc-ccntroids of each of the three Balmcr hncs 
are uncorrelated. For instance, the standard deviation of 
the average velocity as computed from an average of all 
three lines decreases by a factor (ct/3 + ct^ + 0-5)2 /3 com- 
pared to that associated with a single line, as would be 
expected if the errors were mutually independent (here, 
<^f).~i.s arc the respective standard deviations in the ve- 
locities as derived from H/?, H7, and H5). The velocity 
differences between the lines also offer an independent es- 
timate of the measurement uncertainty; e.g. one expects 
the ratio {vfj — v-^) / {a^ + a^) 2 to have a mean of zero and 
a standard deviation of unity if the errors associated with 
each line arc uncorrelated. Our values of the standard de- 
viation of the above ratio for the two nights are 1.2 and 
1.3, i.e., roughly consistent with unity. We deem our error 
estimates to be credible. 

The Fourier Transform of the velocity curve (Fig. 3b) is 
striking in that it shows no strong peak at any frequency, 
not even at the frequencies corresponding to the dominant 
variation in the light curve. We will explore the possible 
cause of this result in later sections. We can, nevertheless, 
place interesting upper limits on the velocity to flux ratio 
of each observed mode. The motivation for doing so is that 
the velocity to flux ratios {Rv) have, to date, only been 
determined for one star (van Kerkwijk et al. 2000); yet, 
these are essential for comparison with theoretical predic- 
tions. While the noise level prevents the velocity curve 
from being used in the absence of external information, 
the light curve provides just such external, independent, 
information as to the periodicities we expect to find in the 
velocity curve. We can exploit this additional information 
by imposing the frequencies wc find in the light curve on 
the velocity curve and measuring the velocity amplitudes 
at these pre-specified frequencies. We can subsequently 
ask if a velocity peak at a known frequency is signficant. 
We detail this procedure below. 

As just described, we looked for modulations in the 
velocity curve by fitting the velocity curve with sinusoids, 
the frequencies of which were fixed at those obtained from 
the light curve. As with the light curve, the calibration rel- 
ative to HS 0507A removed all slow variations and only a 
constant offset was additionally included in the fit. The re- 
sulting velocity amplitudes and phases are listed in Table 
1. We find marginally significant velocity variations at Fl, 
F2, and F6. If these are real, it is surprising that we do not 
find significant velocity variations at F3 and F4, in spite 
of these modes having stronger flux modulations than Fl 
and F6. For F4, this is due in part to the proximity of the 
combination mode F4-I-F6— F5. Excluding all combination 
modes from the fit (theoretically, these are not expected to 
have associated velocity variations) results in an increase 
of the velocity amplitude of F4 to 0.8kms~^ {Ry = 7±7) 
while the amplitudes for the other real modes change by 
less than 0.2kms~^. 



We carried out Monte Carlo simulations in order to as- 
certain the likelihood that the peaks as large as those we 
detected might occur simply by chance. Our simulations 
were conducted thus: we randomly distributed the veloc- 
ities with respect to the times and then fit these velocity 
curves in exactly the same way as the observations i.e., by 
fixing the frequency of the sinusoids to those derived from 
fits to the light curve as described above. This procedure 
was repeated 1000 times and the number of peaks with 
amplitudes larger than those measured at the frequencies 
corresponding to the real modes were counted. The re- 
sults corroborated our error estimates: the modulation at 
F2 was most significant, having a probability of only 2% 
of being a random occurrence, while the modulations at 
Fl and F6 had probabilities of 9% and 4% respectively of 
being chance occurrences. In what follows, we treat these 
measurements as upper limits. 

In summary, wc find that at the frequency of F2, the 
mode with the largest photometric amplitude, there is ev- 
idence for an associated velocity signal with an amplitude 
of 2.6±1.0kms~^, which our Monte-Carlo simulation sug- 
gests has a probability of only 2% of being due to chance. 
If taken at face value, we have thus detected velocity vari- 
ations in a second ZZ Ceti type pulsator. A more con- 
servative conclusion is that the upper limit to the velocity 
variations is 3kms~^ (at 99% confidence). Before address- 
ing what one can infer from the (limits to) the amplitudes, 
we first present a discussion of the chromatic amplitudes, 
from which we attempt to infer the spherical degree of the 
real modes. 

6. Chromatic Amplitudes 

As was noted by Robinson et al. (1995), the wave- 
length dependence of the fractional pulsational amplitude, 
AF\/F\, can potentially provide a clean measurement of 
the spherical degree (£) of the modes as it is not only in- 
dependent of the inclination and m (and therefore also to 
unresolved modes having the same £ but different ni, but 
is also insensitive to the details of flux calibration.^ 

We have computed these fractional pulsation ampli- 
tudes at each wavelength (hence the name, 'chromatic am- 
plitudes') by fitting the real and combination modes listed 
in Table 1 in each 3 A bin. where the choice of bin size re- 
flects a compromise between adequate signal to noise and 
resolution. The frequencies of the modes were fixed at the 
values shown in Table 1 while the amplitudes of the modes 
were left free to vary. At first, we also left the phases free, 
but found that these did not show any significant signal. 

^ Note that the inference of £ from the observed chromatic 
amplitudes only, can, in principle, be made without even the 
qualitative use of models. Although the use of models neces- 
sarily brings with it a dependency on the assumed convective 
efficiency, the qualitative use of model chromatic amplitudes 
computed with differing mixing length prescriptions does not 
affect the mode identification. The £ = 2 mode in Clemens 
et al. (2000) is always better matched with an £ = 2 model 
regardless of the mixing length prescription. 
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F4+F6 27Bs 




1.6 



4000 4500 5000 5500 

A(A) 

Fig. 5. Chromatic amplitudes for the 7 real modes and 2 of 
the stronger combination modes. The period is indicated next 
to the name. From photometry, we know that Fl and F2 are 
modes having the same £ but different m values; likewise for 
F5 and F6. Note the shape and curvature of the continuum 
between the line-cores. The chromatic amplitudes confirm that 
at least the stronger modes share the same £ value (see also Fig. 
7). 



Since Clemens et al. (2000) found from their high signal- 
to-noise ratio data that the phases show very httle vari- 
ation^, we decided to determine the amphtudes with the 



^ Note that Clemens et al. (2000) did find phase changes 
across the Balmer lines, but at levels too low for us to detect. 
They also found a slight, unexplained slope of phase as a func- 
tion of wavelength. We also find a similar variation (2-3°) for 
the strongest modes. 




4000 4500 5000 

A(A) 



5500 



Fig. 6. Comparison between observed and model chromatic 
amplitudes for the strongest mode of HS 0507B (F2, dots) and 
that of ZZ Psc (Fl, full line) shown over a 2000 A range only for 
clarity. Fl of ZZ Psc was identified as having ^ = 1 by Clemens 
et al. (2000). The model chromatic amplitudes overlaid on the 
data are for Toff = 12000 kK, logp = 8 and ML2/a = 0.6. 
All amplitudes have been normalised such that they equal 1 at 
5500 A. The data for ZZ Psc are from Clemens et al. (2000). 



phases fixed to the values listed in Table 1. The resulting 
chromatic amplitudes are shown in Fig. 5 (allowing the 
phases to vary produces a very similar figure). 

Clemens et al. (2000) were able to see a clear contrast 
between the chromatic amplitudes for their £ — 1 modes 
and that for the one i = 2 mode present in their data for 
ZZ Psc. We see no such clear contrasts for HS 0507B (Fig. 
5). Mostly, this simply reflects the fact that HS 0507B is 
about 8 times fainter, and that our signal-to-noise ratios 
are correspondingly lower. In Fig. 6, we compare the ob- 
served chromatic amplitudes for the strongest modes in 
HS 0507B (F2) and ZZ Psc (Fl, £ = 1). We overlay model 
chromatic amplitudes for comparison. These were com- 
puted using model atmospheres (kindly provided by D. 
Koester; earlier versions described in Finley et al. 1997). 
The similarity of the displayed chromatic amplitude of 
HS 0507B and ZZ Psc and the dissimilarity of either to 
the models is striking. Pending improvements to the mod- 
els, we must invent other methods to distinguish between 
£ = 1 and £ = 2 modes. 

A careful inspection of Fig. 5 reveals that F2, F3, F4, 
and F6 look qualitatively similar to each other while the 
continuum shape and curvature - especially between H/3 
and H7 - are somewhat different for F7. Although the 
continuum between the line cores is not smooth for the 
weaker multiplets, the general shape of the line cores is 
not unlike that of the stronger modes. We focus upon these 
differences - the same ones found for ZZ Psc - in order 
to glean any clues on the spherical degree of the pulsation 
modes. 

For ZZ Psc, two distinguishing features of the £ = 1 
modes compared to the one £ = 2 mode were that the 
latter had a larger curvature in the continuum regions be- 
tween the line cores and that the overall slope was steeper. 
In the hope of separating out modes having different £ val- 
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ues, we attempt to measure these two quantities by fitting 
a 2"'^-order polynomial between the lines cores as differ- 
ences in these regions are readily apparent. We use ZZ Psc 
as a test case as at least one t = 2 mode was identified by 
simple inspection. 

In order to minimize the effect of local amplitude vari- 
ations on the curvature and to minimize the covariance 
between the various parameters, we fit y = yo + yiz + y2z'^ , 
where y = In a is the (natural) logarithm of the amplitude 
a and z = In A — (In A) , with (In A) the average of In A in 
the wavelength region of interest. Here, we only show the 
results for the region between H/3 and H7 (4370-4820 A). 
We also define yom = (2/0 - (lna)55oo)/((ln A) - In 5500). 
In this scheme, j/om can be seen to be a measure of the 
slope of the entire spectrum, while yi measures the local 
slope. Plotting 2/2 against yom should group together the 
modes that have similar shapes and curvatures. 

Indeed Fig. 7 shows that F4 from the Clemens et al. 
(2000) data set clearly stands out from the cluster of ^ = 1 
modes. The situation is somewhat less clear-cut for HS 
0507B: while F2, F3, F4, and F5 are consistent with each 
other and with Fl and F2 of ZZ Psc, the position of F6 
is puzzling, although its value of yom is consistent with 
that of the other strong modes. This, together with the 
general appearance of the chromatic amplitudes implies 
that F2, F3, F4, and F5 have ^ = 1 and as Fl and F6 are 
members of known £ = I multiplets, they too must have 
i = 1, even though this is not obvious from their respec- 
tive locations in Fig. 7. The difference in appearance of F7 
mentioned above is borne out by the measured slope and 
curvature. Given its low amplitude, we cannot unfortu- 
nately rule out the contribution of random noise. We note 
with interest that Handler (private communication) finds 
a peak of ~ 0.19% at 3489.09 /iHz i.e. at approximately 
the same splitting as that observed for the 355 s and 445 s 
multiplets suggesting that F7 may be part of an i = 1 
triplet. 

We can check the reliability of the above method by 
repeating the procedure for those combination modes for 
which the £ character of the combination is known. The 
spherical harmonic character of combination modes can 
be deduced from selection rules whenever the i and m 
values of the constituent real modes is known. Thus for 
combinations involving two m = lorm = — 1 modes (for 
£ — 1), the resulting wavelength-dependent pulsation am- 
plitude should have a surface distribution described by a 
component and its position in Fig. 7 should coincide 
with F4 of ZZ Psc (which has £ — 2). Unfortunately, we do 
not have any independent knowledge of m for any of the 
modes in ZZ Psc and cannot use its combination modes 
to verify the locations of the real modes in Fig. 7. Of our 
combination modes, F2-t-F3 is the largest amplitude one 
that satisfies the above m- requirement. Its location in Fig. 
7 is roughly consistent with it having an£ = 2 component, 
although the uncertainties are too large for the measure- 
ment to be significant. In general, our conclusion from 
the above exercise is that although one obtains quantita- 
tive measures for the differences between modes having 
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Fig. 7. Measures of the observed (HS 0507, ZZ Psc) and model 
curvatures in the chromatic amplitudes between H/3 and H7. 
Energy transport by convection is comparable in the MLl / a— 2 
and ML2/q=0.6 versions. The model amplitudes were con- 
volved with a Gaussian having a FWHM of 4.6 A to match 
the seeing profile and extended over the 8. 7 slit. 



differing £ indices, it merely serves to corroborate the con- 
clusions arrived at by visual inspection of the chromatic 
amplitudes. 

As an aside, we note that using all the real modes listed 
in Table Al in van Kerkwijk et al. (2000) results in one 
other mode, F8 (920 ± 3 s), being almost coincident with 
the position of the £ = 2 mode (F4, 776 s) in Fig. 7, while 
FA (500 s) is consistent with the position of F4 within 
the errors. Taken at face value, this implies that the 920 s 
mode also has £ = 2. Kleinman et al. (1998) assumed £ — 1 
for their 918 s mode so it would be interesting to check 
the validity of this assumption by assigning a value of 
£~2to this mode in the pulsation models and attempting 
to quantify the resulting changes, if any, to the derived 
properties of ZZ Psc. Given the relatively low amplitude 
of F8 (0.47%), other independent constraints would, of 
course, be desirable. 

We repeated the above procedure on the model chro- 
matic amplitudes for a range of temperatures and a num- 
ber of different sets of mixing- length parameters. We show 
two examples in Fig. 7 one for ML2/a = 0.6, which was 
found to yield the best description of the average optical- 
to-ultraviolet spectra (Bergeron et al. 1995) and one for 
MLl/a = 2.0, which lies in the parameter space of models 
with which the observed magnitude difference between HS 
0507A and B was best reproduced (Jordan et al. 1998); see 
Jordan et al. (1998) for a description of the terminology. 
We find that all models, like the two shown, have sys- 
tematically higher curvature than the observations. Thus, 
we extend to all models the conclusion of Clemens et al. 
(2000) that while the salient features of the chromatic am- 
plitudes are reproduced by the models, the details leave 
much to be desired. It may well be that better agreement 
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will only be reached once a better description of convec- 
tion becomes available. 

7. Comparison with Pulsation Theory 

In this section, we attempt to place our observations 
within a theoretical framework by comparing various 
quantities derived from our observations to those derived 
from theoretical considerations by Wu & Goldreich (1999), 
Goldreich & Wu (1999a), and Wu (2001). Our primary aim 
in this section is to check whether theoretical expectations 
arc borne out by the observations for both real and com- 
bination modes. Our secondary aim is to use our observa- 
tions to perform consistency checks on some of the theo- 
retical parameters used in the above studies. Specifically, 
these parameters are the thermal time constant of the 
convection zone (tc^), and the parameter (|2/3 + 7I) that 
depends on the depth of the convection zone. We addi- 
tionally check whether the inclination angle derived by 
Handler fc Romero- Colmenero (2000) is consistent with 
our data. Before proceeding, we stress that any compar- 
ison with the theory is rendered difficult by a number of 
factors, in particular, by unresolved multiplets and by the 
presence of real and combination modes of low amplitude. 
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8. The relative velocity to flux amplitudes (a) and pliases 
(b) for HS 0507B and ZZ Psc. The values for ZZ Psc are taken 
from Table 1 in van Kerkwijk et al. (2000). 



7.1. Real Modes 

Following van Kerkwijk et al. (2000), we measure the rel- 
ative fliix and velocity amplitudes and phases of the real 
modes with Ry = Ay/(27r/^L) and A$y = $y - 
respectively. In Table 1, we list the values of Ry for all 
real modes. Here, one should bear in mind that the detec- 
tion of even the strongest mode is only marginal. For the 
weaker modes, the Ry values represent upper limits. 

As a ZZ Ceti type white dwarf cools, the depth of 
the convection zone increases and longer period modes 
are excited. However, the flux perturbations are also more 
effectively attenuated so the emergent flux variations at 
the photosphere for a mode of fixed internal pulsation 
amplitude are reduced. Now turbulent viscosity in the 
convection zone ensures negligible vertical velocity gradi- 
ents, making the horizontal velocities effectively indepen- 
dent of depth within the convection zone (Brickhill 1990; 
Goldreich & Wu 1999b). Thus at a fixed frequency, Ry 
should be smaller in white dwarfs with thinner c;onvec- 
tion zones, i.e. the hotter pulsators. Similarly, for differ- 
ent modes within the same object, the flux attenuation 
increases with increasing mode frequency, while velocity 
variations pass undiminished through the convection zone. 
Thus Ry is expected to increase with increasing mode fre- 
quency. A^y is expected to be equal to 90° for the adia- 
batic case (flux leads maximum positive velocity by tt/2) 
and to tend towards 0° with increasing frequency, as the 
flux variations are increasingly delayed. 

In Fig. 8, we show the observed values of Ry and 
A^v for both HS 0507B and ZZ Psc. The observations 
are broadly consistent with expectations, in that the av- 
erage values of Ry are reasonable, and values of A$y are 



between and 90°. We find no clear evidence, however, 
for the expected trends with frequency as the uncertainties 
are too large. 

The longest period observed mode can be used to de- 
duce a rough lower limit to the value of Tcq, the thermal 
time constant of the convection zone, as the overstability 
criterion loTc„ ^ 1 must be satisfied for driving to exceed 
damping (Goldreich & Wu 1999a). Using F4 from Table 1 
yields a value for of ^ 118 s. 

7.2. Combination Frequencies 

In models, the sharp max;ima and shallower minima that 
typify the light curve are the result of the variation in 
the depth of the convection zone in response to the per- 
turbations; this variation distorts the flux variation and 
gives rise to linear combinations of frequencies when trans- 
formed into frequency space (Brickhill 1992; Wu 2001). 
Thus, combination modes are not expected to have asso- 
ciated physical motion. 

The relative photometric amplitudes and phases of the 
combination modes with respect to their constituent real 
modes can be expressed as Rc = A^^-' / {rnjA^A'j^) where 
Tiij is the factorial of the number of different modes in 
the combination for 2-mode combinations and A^c = 
<I>^-' — (<I>|^±$^) respectively. Relations linking theoretical 
quantities to observed ones are provided by Eq. (15) and 
(20) in Wu (2001). We quote these here for convenience: 

A$c = tan-i f . ,\ ) (1) 
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|2/3 + 7|(u^»±o^j-)rco ^Ui, 

Here, w is the angular frequency of the (real) mode, /? 
and 7 are related to the entropy variation, the response 
of the superadiabatic layer to the pulsation, and to the 
thermal relaxation time at a given depth in the convection 
zone. G^'^.^^g^'g™" is the geometrical multiplication 
factor which depends on the inclination, 6o between the 
pulsation axis and the line of sight^ and ay ~ 0.4 is a 
factor relating photometric amplitudes in the V band to 
bolometric ones. 

In what follows, we take the inclination angle between 
the pulsation axis and the line-of-sight to be 79° as in- 
ferred by Handler & Romero- Colmenero (2000). We cau- 
tion, though, that its derivation assumed not only that 
the pulsation axis is closely aligned to the rotation axis 
of the star, but also that multiplets have the same intrin- 
sic amplitude. While the former assumption is probably 
justified, the latter is not self-evident as the amplitude of 
the modes may depend highly non-monotonically on fre- 
quency (especially if they are determined by parametric 
instability e.g. Wu & Goldreich 2001). Observationally, if 
the intrinsic amplitudes of the multiplet components are 
identical, one would expect the m = ±1 components of the 
triplets to have the same amplitude, and the ratio of the 
m = to the m = ±1 to be the same for different multi- 
plets. While the observations are generally consistent with 
the above, there is a glaring exception: the m = ±1 com- 
ponents of the 355 s triplet, which includes the strongest 
mode, have very different amplitudes both in our data 
(Fl, F2), and in those of Handler & Romero-Colmenero 
(2000). The derived value of ©o should therefore only be 
considered as approximate. 

Figure 9 shows the observed and theoretically expected 
values of Rc- For the latter, we used 0o = 79°, Tcq = 118 s 
(as derived from the longest period ovcrstable mode) and 
1 2/3 -|- 7 1 = 12.6, as estimated from theoretical arguments 
by Wu (2001). Given the numerous assumptions that have 
to be made in deriving theoretical values of Rc and the 
uncertainty in the measured mode amplitudes, especially 
at lower frequencies, the agreement is adequate: values of 
Rc arc expected to be lower for lower frequency combina- 
tions (wTco — 1) and almost constant for higher frequency 
combinations (wTcq 1) as can be seen from Eq. (2). The 
best agreement with expected values of Rc is shown by 
combinations involving two m = — 1 modes, these are also 
the largest amplitude modes in general. 

For combinations having rrii irrij = ±2; ±1, the theo- 
retical value of Rc is independent of the inclination of the 
pulsation axis to the line-of-sight (8o) and depends only 
on the assumed value of 1 2/3-1-7 1, while for m, = —nij = ±1 

^ Values of this factor for grey atmospheres and £i = 
tj = I are listed in Table 3 of Wu (2001). It is equal 
to 0.65 -I- O.45/cos^0o for rrn = mj = combinations, 
0.65 + 0.90/ sin^ &o for rrii = —rrij = ±1 combinations, and 
0.65 for all others. 



U.U. u. u. b.1 1 1 1 1 1 U.I I 1 1 b. U. 

Ill + r<i++++ +++ + + 




Frequency i^lHz) 

Fig. 9. Measure of strength of two-mode combinations, Rc ~ 
A^j^^ /{riijA'iA'j^), as a function of frequency. Filled and open 
squares denote those combinations for which values of m are 
known (from Handler & Romero-Colmenero 2000) and un- 
known respectively. Theoretically expected values of Rc for 
£i = (.j = 1 are calculated using Eq. (2) with an inclination 
angle of 79° and Tcq = 118 s. The assumed value for |2/3-|-7| is 
indicated. The upper and lower arrows denote the asymptotic 
value of rrii = —rrij = ±1 combinations with Go = 50° in Eq. 
(2), with |2/3 -I- 7I = 12.6 and 9 respectively. It is an example 
of how other values of Go would shift the theoretical curves. 
(As indicated in the text, the mi ± rrij = ±2, ±1 curve is in- 
dependent of Go). The dashed line ending in the open triangle 
shows where the F1-I-F5 combination would be shifted to if it 
has been misidentified (see text and Table 1). 

combinations, Rc is dependent on |2/3 -|- 7I and ©0, but 
is independent of for higher frequency combinations 
(wTco » 1). Taken together, /3 and 7 represent a mea- 
sure of the depth of the convection zone as a function 
of effective temperature of the white dwarf and thus re- 
late directly to the width of the instability strip: the lower 
|2/3-|-7| is, the wider the instability strip (Wu 2001). Figure 
9 shows that better agreement with the observed values 
of Rc is obtained using a somewhat lower value of 9 for 
|2/3 -|- j\. For this value, the rrii = — = ±1 modes 
are also better reproduced. This supports a high value 
for the inclination angle: for smaller ©0, the theoretically 
expected Rc values for rui = —mj = ±1 would shift up- 
wards and away from the observations. 

The relative phases are shown in Fig. 10. These are ex- 
pected to be in phase with their corresponding real modes 
i.e. A'I'c' ^ 0. While many combinations follow the ex- 
pected trend, there are some large deviations; we discuss 
this further below. 

In the above, we used Tc„ ~ 118 s, as derived from the 
longest period overstable mode. In principle, one might 
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Fig. 10. Observed (filled squares) values of A$c and the theo- 
retically expected values as a function of frequency, as obtained 
from Eq. (1), assuming a value of 118s for Tcq. As in Fig. 9, 
the dashed line ending in an open triangle, shows the position 
of F2+F6. 

improve on this by fitting the whole set of combination 
modes for Tcq, \2(3 + ^\ and 6o simultaneously. However, 
Wu (2001) showed that for pairs of sum and difference 
combination frequencies arising from the sam,e real modes 
(e.g. F1±F5), the ratio of the amplitudes does not depend 
on |2/J+7| and thus one could derive directly. We quote 
Eq. (21) from the above-cited paper: 

(c..+o.,)vA + (^,-^,)^T;^gr-r^ 

Note that the geometrical factor and hence the depen- 
dence on 6o cancels out if at least one constituent real 
mode has m = 0. 

From the list of combination frequencies in Table 1, 
there are three pairs that may be used, viz., F1±F5, 
F3±F4, and F2±F4. None of these have an to = com- 
ponent, but we can still solve Eq. (2) for r,,-, using the 
estimate of Oq = 79° of Handler & Romero-Colmenero 
(2000). The results are listed in Table 2; we see that the 
values of Tc„ thus derived are compatible with the value 
(~ 118 s) derived from the longest period real mode ex- 
cited. Using 6o = 79° together with Eq. (2) and 3 yields 
a value of Tc„ that is entirely consistent with the value 
found from the longest period real mode. This agreement 
lends some confidence to Eqs. (2) and (3). We therefore 
point out that only the choice of m = — 1 for F4 resulted in 
physically acceptable solutions; this is consistent with all 
observed multiplet components in our data having to 7^ 
and with the indication from the chromatic amplitudes 
(Sect. 6) that all modes that give rise to combinations 



have ^ = 1. As we point out below, the F1-I-F5 combina- 
tion may in fact be F2-I-F6. In view of this, it is hardly 
surprising that the F1±F5 pair came to naught. 

Table 2. Values of Tcq derived from combination pairs 



Combination pair 


Too (s) 




F1±F5 






F2±F4 


112 {rn = 


-1) 


F3±F4 


122 (rn = 


-1) 



Notes. For the F2,3±F4 pairs, £j = l,m = 1,0 resulted in 
unphysical solutions. Likewise for the F1±F5 pair. All possible 
permutations of £i = l,mi = = 2,mj = ±2, ±1,0 gave 

unphysical solutions for the F2±F4 pair. Likewise for F3±F4 
except for ij = 2,mj = —1, for which Tcq=19s. 

In Figs. 9 and 10, there are a number of combinations 
for which Rc and A^c are rather discrepant from the the- 
oretically expected values. These discrepancies might sim- 
ply be due to unresolved, low amplitude modes, to which 
the measured phases in particular are rather sensitive. 
However, some of the anomalies can be explained more 
easily, as (i) resulting from degenerate combinations^ that 
are an inevitable consequence of equal frequency split- 
tings, and (ii) the possibility for low amplitude real modes 
to have high amplitude combinations. We discuss each of 
these possibilities below, using our observed combination 
modes as examples. 

A case in point is our F1-I-F5 combination that lies 
very close in frequency to the F2+F6 combination. In 
hindsight, the anomalous value of Rc points to a misiden- 
tifi-cation, but only because we now have an £ value for 
the Fl and F5 modes, bolstered by an estimate of the in- 
clination angle. Including F2+F6 in the fit of the Fourier 
Transform of the light curve (Sect. 4) instead of F1-I-F5 
results in a value of Rc (8 ± 1) that is more compatible 
not only with the theoretically expected value, but one 
that is also more in line with other combinations involv- 
ing TOj = l,mj = —1 (see Fig. 9). The value of A^c, 
—28 ±5°, is also in better agreement (sec Fig. 10).^ 

More problematic are combinations that are degener- 
ate with TOj = TOj = combinations, as with ©0 = 79° 
these combinations have very high Rc values (sec Fig. 9), 
as a result of strong cancellation for the real mode and lit- 
tle cancellation for the combination (which has an ^ = 
component). This means that even real modes with weak 
observed amplitudes can give rise to relative strong com- 
binations^'^. As an example, consider our F1-I-F3 {rrii = 

* Degenerate combinations are a particular problem for HS 
0507B, since the frequencies of several real modes have near- 
integer ratios e.g. F3:F4 ~ 4:3 (see Table 1 and Handler & 
Romero-Colmenoro 2000) . 

^ Our justification for initially preferring F1-I-F5 is as stated 
in Sect. 4 i.e. that we obtained a slightly larger fitted amplitude 
and a correspondingly lower x?ed foi^ F1-I-F5. 

Indeed, for 60 = 90°, the combination mode would be 
visible oven though the corresponding real modes would not! 
Obviously, this could lead to gross misidentifications. 
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l,mj = —1). This combination has the same frequency 
as f3+f9 i.e. 1796.70 + 2814.02 /uHz, nn = rrij = - of 
Handler & Romero- Colmenero (2000). Fitting our Fourier 
Transform using the frequencies of Handler & Romero- 
Colmenero (2000) (Sect. 4), we find that the sums of the 
phases of the two combinations of real modes arc almost 
the same: -29° for f3-|-f9 and -32° for F1-I-F3. Hence, 
these combinations will add coherently. If the real modes 
in the different multiplets have the same intrinsic ampli- 
tude, one would expect both combinations to have roughly 
the same observed amplitude, and hence this could plau- 
sibly lead to an anomalously high value of Rc for F1+F3. 

Another such combination is F4:m=-i+F6m=i 
(1346. 1-1-2248. 4^Hz) identified by Handler (private com- 
munication) as /2„=_i + /4„=i (1793.29-hl800.7/LtHz). 
It has an amplitude that is not only the largest among 
the combination modes but is also larger than that of 
2 of our real modes (Fl and F7). Another combination 
that also adds up to ~3594/xHz is the first harmonic 
of the m = — 1 component of the 556 s multiplet i.e. 
2 X 1796.7 or 2/3 in Handler & Romero- Colmenero 
(2000). As mentioned earlier, even though the m = 
components are weak, their harmonics are not. This, 
coupled with the integer frequency ratios for the modes 
involved means that different combinations can (and do) 
have degenerate frequencies. In this respect, it may be 
interesting that the one mode from the study of Handler 
& Romcro-Colmcncro (2000) that changed in amplitude 
is /4„^i (Fig. 4): perhaps this is the result of a resonance 
between these modes. This could explain the anomalously 
high value of Rc for the F4-I-F6 combination. 

In summary, therefore, we find that most of the appar- 
ently discrepant values of Rc and A$c can be explained 
by the contribution of combinations of unresolved or low- 
amplitude modes, degenerate combinations, and low fre- 
quency noise. 

8. Squaring the circle 

We have used high signal-to- noise, time-resolved spectra 
in an attempt to identify the spherical degree of the pul- 
sation modes and to infer some of the properties of the 
outer layers of a pulsating white dwarf. In the process of 
doing so, we have tested several aspects of the very theory 
that was used in the analysis. Real and combination modes 
were used in conjunction with wavelength-dependent frac- 
tional amplitudes and velocity amplitudes in order to piece 
together a consistent picture. At almost every step, our re- 
sults were compared with those of ZZ Psc, a brighter and 
better-studied white dwarf. 

Our measured frequencies and amplitudes (sec Table 
1) for all modes except one are in agreement with those 
reported by Handler & Romero- Colmenero (2000). 

We measured the linc-of-sight velocities associated 
with the pulsations and found a marginally significant 
modulation at the frequency of the strongest mode. This 
can be taken as a detection of surface motion in a sec- 
ond ZZ Ceti type pulsator, or, more conservatively, as a 



stringent upper limit to such motion. From all modes, we 
find that the average value of (or upper limit to) the ratio 
of the velocity to flux amplitudes, i?y, is just over half 
of that observed in ZZ Psc. The difference in the thermal 
time constant of the convection zone, Tcq , for the two white 
dwarfs of about a factor of two (118 s c.f. 250 s, Wu 2001) 
translates into a factor of two in Ry as the latter roughly 
scales with (Goldreich & Wu 1999a,b). Our measured 
values of Ry are therefore consistent with expectations. 

Based on the general appearance and shape of the 
continuum between the line cores of the chromatic ampli- 
tudes, we deduced that modes Fl - F6 were consistent with 
having £ = I. This adds some confidence to the pulsation 
models of Kleinman et al. (1998) that assign £ = 1 iden- 
tifications to modes having periods of 284 s, 355 s, 445 s, 
and 555 s. 

The chromatic amplitudes of these £ = 1 modes were 
found to be more similar to the £ = 1 modes in ZZ Psc than 
to those expected from the models, supporting the conclu- 
sions of Clemens et al. (2000). In order to make our results 
more quantitative, we devised a measure of the curvature 
and slope of the chromatic amplitudes and used this in an 
attempt to separate £ = 1 and £ — 2 modes for both the 
observations and the models (Fig. 7). Although this proce- 
dure did not alter our mode idcntifcation based on simple 
inspection for HS 0507B, it did serve to quantify differ- 
ences between the models and observations. Furthermore, 
we found that it worked very well for ZZ Psc and that it 
yielded two additional potential £ = 2 modes, at 920 and 
500 s. 

The relative amplitudes of the combination modes 

proved to be a useful tool, as these in general follow ex- 
pected trends for given £ and m. Exploiting this, we at- 
tempted to redetermine using combination frequency 
pairs, and obtained good agreement with the value de- 
rived using the longest period real mode. A by-product of 
this exercise was an indication that F4 has m = — 1 as no 
other m value - including those combinations arising from 
£i = l,£j =2 combinations ~ yielded physically accept- 
able values of , thereby supporting our £ identification 
for F4. Thus, combination frequencies can potentially pro- 
vide indirect constraints on £ and m values even when no 
splittings are observed. Taking at face- value, we infer 
that HS 0507B has a shallower convection zone than ZZ 
Psc, consistent with HS 0507B being slightly hotter. 

The phases of the real and combination modes were 
slightly more difficult to reconcile with theoretical expec- 
tations. For combination modes this is likely to be due 
to unresolved modes that plague a relatively short time 
scries. Data from longer, uninterrupted time scries, for in- 
stance, from the Whole Earth Telescope is probably better 
suited to such analysis. We showed that weak, unresolved 
modes are probably being manifested in the amplitudes 
and phases of the combination modes. For the real modes, 
although our measured values of A$y are less than 90°, 
the trend with frequency is opposite to that expected; this 
is also true for ZZ Psc (van Kerkwijk et al. 2000). 
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Potential future work could involve focusing on using 
the observed chromatic amplitudes to calibrate the tem- 
perature stratification in current white dwarf atmosphere 
models and may provide insight into areas where these 
models fail to reproduce the observations, while attempts 
to match the stringent constraints of the observed period 
structure and ^ and m identifications using pulsation mod- 
els may yield a unique asteroseismological solution for HS 
0507B. 
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Appendix A: Corrections for instrumental flexure 
and differential refraction 

We were especially fastidious in correcting the wavelength 
scale for the effects of flexure and differential atmospheric 
refraction. We measured the instrumental flexure from the 
positions of the Oi A5577A sky line which were derived 

by cross-correlating the flux-calibrated spectra using an 
average of all spectra as a template. We found the posi- 
tions to be adequately represented by a third-order poly- 
nomial fit (Fig. A. la). The effect of differential refraction 
on line and continuum intensities is non-negligible even 
at modest airmasses and increases rapidly with decreas- 
ing wavelength. Having corrected for instrumental flexure, 
the resulting positions of the Balmer lines vary linearly 
with sin a tan Q where a is the difference between the par- 
allactic and position angles and C, is the zenith distance. 
This is indicative of the remaining shifts being dominated 
by refraction. (Upper curve Fig. A. lb). We computed a 
value for the refractive index relative to that at H7 us- 
ing the standard prescription of Stone (1996) and use this 
value to correct all pixels in the spectrum. As the position 
of the stars in the slit was altered after a correction for 
a guider malfunction (on the second night) compared to 
that at the beginning of the observing run, we applied a 
5.1 A shift to the spectra of both components. This shift 
was computed by measuring the wavelength of H7 of HS 
0507A in a narrow-slit (0'.'7) frame and insisting that the 
average wavelength of H7 as obtained from line-proflle 
fltting - be the same. For complete consistency, we applied 
a 1.74 A shift to the spectra of the previous night. 

Appendix B: HS 0507A 

As was mentioned in Sect. 3, the non- variability of HS 
0507A is crucial to the velocity measurements and the 
continuum variations of HS 0507B as these are carried 
out with respect to HS 0507A. The measurements of HS 
0507A are, of course, not carried out relative to any other 
star, this is further elaborated in Fig. B.l. 
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Fig. A.l. (a) The lowest curve (offset by 3.5 pixels) shows 
the pixel positions of the Oi A 5577 A line as a function of 
time. The upper curve (1) shows the variations in the position 
of H,5 while (2) shows the variations in H/3 after correcting for 
instrumental flexure, (b) The dots show the resulting positions 
of H/3 after shifting with respect to H7 and the pluses the final 
positions of H/3 after correcting for both flexure and refraction 
(offset by 5 pixels), (a) and (b) together show the various stages 
in the processing that lead to the flnal positions of the lines. 
Note that although the above figure is for the first night only, 
it is representative of the procedure for both nights. 
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Fig. B.l. (a) Fourier transforms of the light curves of HS 
0507B (top) and HS 0507A (offset by -0.3%) (b) Fourier 
transforms of the velocity curves of HS 0507B (top) and HS 
0507A (offset by -3kms"^). The average velocity curve of HS 
0507A, as used here, was constructed from the flexure and 
refraction-corrected positions of the three strongest Balmer 
lines. Polynomial fits were used to describe wandering in the 
slit for HS 0507A whereas H7 of HS 0507A was used to cor- 
rect for wander of HS 0507B as described in Sect. 3. The er- 
ror in a single measurement of the velocities for HS 0507A is 
■~ 9kms~^. Note the lack of significant velocity signal in HS 
0507A at the frequency of the strongest mode as obtained from 
the flux variations of HS 0507B. The Fourier transforms shown 
here are for the first night only. 



